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We analyze by means of extensive computer simulations the out of equilibrium dynamics of 
Edwards- Anderson spin glasses in d = 4 and d = 6 dimensions with ± J interactions. In particular, 
we focus our analysis on the scaling properties of the two-time autocorrelation function in a range of 
temperatures from T = 0.07T C to T = 0.75T C in both systems. We observe that the aging dynamics 
of the ±J models is different from that observed in the corresponding Gaussian models. In both 
the Ad and 6d models at very low temperatures we study the effects of discretization of energy 
levels. Strong interrupted aging behaviors are found. We argue that this is because in the times 
accessible to our simulations the systems are only able to probe activated dynamics through the 
lowest discrete energy levels and remain trapped around nearly flat regions of the energy landscape. 
For temperatures T > 0.5T C in 4d we find logarithmic scalings that are compatible with dynamical 
ultrametricity, while in 6d the relaxation can also be described by super-aging scalings. 



I. INTRODUCTION 

After more than twenty years of extensive research, 
the physics of spin glasses is still far from being com- 
pletely understood. The inherent complexity of the phys- 
ical scenario together with unsurmountable mathemati- 
cal difficulties undermine even the simplest theoretical 
approaches that attempt to include basic realistic ingre- 
dients. Notwithstanding that a large amount of informa- 
tion could be extracted from the analysis of the mean 
field model (i.e. the Sherrington-Kirkpatrick model 1 ' 2 ), 
the extension of those results to finite dimensional spin 
glasses remains a matter of hard controversy in the statis- 
tical physics community 3 ' 4 . In this context large scale nu- 
merical simulations emerged as a valuable aid in gaining 
physical insight into more realistic models 5 ' 6 ; however, 
they pose such a serious demand over existent computa- 
tional capabilities that the study of the low temperature 
phases of these systems may be considered in many re- 
spects to be still in an exploratory state. 

Regarding the out of equilibrium behavior of finite di- 
mensional spin glasses, much attention has been devoted 
to models with continuous couplings distributions 6 ' 7 . 
One important feature of these systems is that the ground 
state is unique. For Gaussian couplings, the aging be- 
havior of models in 3d and 4<i reported in the literature 
seems to be compatible with the simple aging or weak 
interrupted aging scenarios (to be defined below). These 
particular scenarios suggest a rather simple phase space 
structure with a unique relevant time scale, namely the 
age of the system. This apparently simple behavior could 
be a consequence of the existence of strongly separated 
temporal scales. Because of this during the time span of 
a simulation or experiment at very low temperatures the 
system could not be able to cross over between different 
time regimes and only one relevant time scale is probed. 



One may wonder whether the relatively simple dynam- 
ical behavior observed in spin glasses with continuous dis- 
tributions at very low temperatures can also be expected 
in systems with discrete couplings, which present a strong 
degeneracy of their ground states as well as a noticeable 
discretization of their low-energy spectrum. That might 
lead to very low temperature effects that could not be 
observable in continuous spin glasses. 

It is known that at temperatures relatively close to T c 
there is no qualitative difference in the phenomenological 
characterization of aging 10 between the continuous and 
discrete couplings systems. However, at lower temper- 
atures, effects deriving from the discrete nature of the 
energy spectrum as well as from the high degeneracy of 
the ground states may become apparent in discrete mod- 
els. 

Scaling properties contribute to a quantitative descrip- 
tion of complex phenomena, even in cases where a general 
theory is lacking 3 ' 11 ' 12 . In this respect, different dynami- 
cal universality classes may emerge for which similar scal- 
ing functions describe the dynamics of different systems. 
The knowledge of these scaling rules gives considerable 
insight into the nature of the underlying dynamical pro- 
cesses. 

In this paper we present the results of an extensive 
numerical study of the aging dynamics and scaling prop- 
erties of the two-time autocorrelation functions for the 
± J Edwards- Anderson spin glasses in dimensions d = 4 
and d — 6 at temperatures which cover the whole low 
temperature phases of the models. 

Our results show important differences between the ag- 
ing dynamics of the ±J models compared with the corre- 
sponding Gaussian models. At the very low temperatures 
T = 0.07T C and T = 0.15T C we find interrupted aging sce- 
narios both in id and 6d. In particular, for the lowest 
temperature we find a very small exponent indicative of a 
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rapid relaxation to a stationary dynamics. At first sight 
this seems hard to reconcile with the expected very slow 
relaxation which takes place at these extremely low tem- 
peratures. Nevertheless this can be seen as a consequence 
of the discrete nature of energy levels which can produce 
some non trivial time dependent phenomena at very low 
temperatures at time scales roughly independent of the 
waiting time t w . The systems simply do not have enough 
time to relax over barriers dependent on t w . So this dy- 
namics is completely different from the usual long time 
non-equilibrium relaxation but is not trivial due to the 
discreteness of the low lying energy levels. At T = 0.5T C 
we observe some differences in the behavior in 4d and 
in 6d. In 4c? logarithmic scalings, which are compatible 
with dynamic ultrametricity, work very well even at this 
relatively high temperature. In 6d the scenario is less 
clear since we observe that although logarithmic scalings 
are very good it is also possible to scale the data with a 
super-aging scaling form. As d = 6 is the upper critical 
dimension the relaxation dynamics should present some 
features of the infinite range model. We find insted that 
the behavior is more similar to that at lower dimensions 
than to the mean field model. 

The paper is organized as follows: in section II we 
define the model studied and the observables measured; 
in section III and IV we analyze the results of the 4d 
and 6d models respectively. Finally the conclusions and 
a discussion are presented in section V. 



II. MODEL AND METHOD 

The system consists of a d-dimensional hypercubic lat- 
tice of Ising spins which interact according to the follow- 
ing Hamiltonian: 



H = 



Jij^iSj 



(1) 



where the symbol indicates that only first neighbor 
pairs i,j are taken into account. The coupling constants 
Jij are binary random variables chosen from the following 
probability distribution: 



p(J ij ) = l(6(J ij -l) + 5(J ij + l)) 



(2) 



The time evolution of the model is governed by a stan- 
dard heat-bath Monte Carlo process with sequential ran- 
dom update. The imposed boundary conditions were pe- 
riodic in 4d and helical in 6d 13 . In the practical imple- 
mentation of the numerical algorithm a significant in- 
crease in speed was accomplished by using multi-spin 
coding 13 so each spin and coupling constant demand just 
one bit of information each for storage. In this way we 
can run many replicas at the same time at the cost of a 
single realization. In all cases the dynamics is initiated 
from a random configuration, simulating a sudden quench 
from infinite temperature into the spin glass phase. 



One straightforward way to characterize the out of 
equilibrium dynamics of complex magnetic systems is 
through the analysis of the two-time autocorrelation 
function C(t,t'), which can exhibit history dependent 
features usually referred to as aging. A system that has 
attained thermodynamic equilibrium will show a station- 
ary dynamics for which only time differences make phys- 
ical sense, and therefore C(t,t') = C(t — t'). However, 
complex magnetic systems such as spin glasses show a 
much more complex behavior due to the presence of an 
extremely slow relaxational dynamics. These materials 
may be out of equilibrium for spans of time longer than 
any available time scale in the laboratory. In this circum- 
stances insight into the ongoing processes can be obtained 
by studying the scaling properties of dynamical quanti- 
ties like C(t w +t,t w ), where the waiting time, t w , stands 
for the age of the system measured after a quench into 
the spin glass phase, and t stands for the time measured 
since the age t w . 

In the numerical experiments we compute the quantity 



C{tnj ~f~ t 7 t w ^ 



1 N 



• t)Si(t w ) 



=1 



(3) 



where we have denoted by [. . .] av an average taken over 
several realizations of the random couplings and thermal 
histories. An additive form for the autocorrelations was 
assumed: 



C(t t 



+ t,t w ) — C s t(t) + C ag (^— -r 



t) 



h(t w ) 



) • (4) 



The stationary part C s t(t) is well described by an alge- 
braic decay of the form C s t(t) = At~ x ( T ^ + q. There is 
no theoretical basis for determining the scaling function 
h(z) appearing in the aging part of the autocorrelations. 
There have been proposed a couple of simple forms which 
work reasonably well for relaxations in the vicinity of sim- 
ple aging, in the sense that the correct scaling variable is 
of the form t/ti^ with /i ~ l 7 . In our case it was not pos- 
sible to find a unique realization of the function h(z) for 
the whole range of temperatures studied. In the scaling 
analysis we restricted ourselves to find the relevant scal- 
ing variables for the aging part of the autocorrelations. 



III. D=4 

The simulations of the four dimensional Edwards- 
Anderson model were done for systems of linear size 
L = 12 (T = 0.07T C ) and L = 10 imposing periodic 
boundary conditions. We recall that the critical temper- 
ature for this model has been estimated to be T c w 2 8 ' 9 . 



A. Very low temperatures: probing the 
discreteness of the lowest energy levels 

In Figure 1 we show the behavior of the two-time au- 
tocorrelation function C(t w + t,t w ) at the temperature 
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T = 0.07T C . The waiting times are t w = 2 fc , k = 14 . . . 19, 
and the simulation was run up to t = 10 7 MCS. There 
is a weak dependence on the wainting time and the sys- 
tem relaxes to a stationary regime presenting a strong 
sub-aging behavior. 
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FIG. 2: Interrupted aging scaling for the asymptotic time 
regime of the data shown in figure 1. The stationary decay 
of the correlations has been subtracted, fit parameters are 
A = 0.04 and x = 0.57. The interrupted aging exponent is 

n = o.i. 



FIG. 1: Two-time autocorrelation function for the ±J spin 
glass in d — 4 with linear size L = 12 after a quench to 
T — 0.07T C in a double logarithmic plot. The waiting times 
are t w = 2 k ,k — 14 . . . 19 from bottom to top. Inset: the 
result of a microcanonical run after t w for the three longest 
tw' s. 

In Figure 2 we show the best scaling obtained for the 
long time behavior of the autocorelation function after 
subtraction of the stationary part: 

C st (T) = At-*W+q. (5) 

In this particular case, the best fit was obtained for 
x(T) = 0.57 (a rather high value) and q w 0.93. In the 
(weakly) t w dependent region it is possible to distinguish 
two regimes: a first one when the system falls out of 
the quasi-equilibrium and a second regime at long times. 
Both are very well described by different interrupted or 
sub-aging scalings of the form: 

C ag (t w +t,t w ) = c(-jr') . (6) 

For very long times, the best fit yields an exponent 
/i ~ 0.1 which characterizes a very strong sub-aging 
regime. It is also possible to collapse the data in the 
time window corresponding to the first out of equilib- 
rium regime (i ~ t w ) with the form given by Eq.(6) with 
an exponent /z ~ 0.27 as can be seen in figure 3. Note 
that the first sub-aging regime disappears gradually as 
t w grows and effectively, after the longest waiting time 
studied t w ~ 5 • 10 5 , these slow modes have practically 
equilibrated. After these very long waiting times the sys- 
tem evolves quickly to a stationary situation with a very 
short characteristic aging time — t w A . The observation of 
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FIG. 3: Interrupted aging scaling for the intermediate time 
regime of the data shown in figure 1. This time regime is char- 
acterized by an interrupted aging exponent /i = 0.27 higher 
than the corresponding one seen in the asymptotic regime of 
figure 2. 



these sub-aging regimes at this very low temperature is 
consequence of the discrete nature of the low lying energy 
levels. In fact the smallest time scale for activation in this 
energy landscape is of the order of exp (2/T) w 1, 6 ■ 10 6 
which is approximately the time where the correlation 
leaves the quasi-equilibrium regime (see Figure 1). Up 
to this time scale the system relaxes effectively in the 
flat energy landscape defined by the sites with zero local 
field which are free to flip (besides a smaller group which 
can still contribute to lowering the energy). Note that 
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the value of the correlation in the plateau C = q in this 
regime does not correspond to the equilibrim order pa- 
rameter qEA but instead to one minus the fraction of sites 
with zero field 21 . After a time of the order of exp (1/0.07) 
thermal activation begins to take place through the low- 
est lying barriers and the correlations decay to zero as 
the system expands its available phase space. As a fur- 
ther test for this interpretation we have performed a set 
of simulations in a microcanonical ensemble which are 
shown in the inset of Figure 1: for the three longest 
waiting times we have allowed the system to relax up 
to t w and from then on only flips of spins with zero local 
fields were allowed. In these conditions the correlations 
decayed only to the plateau suggesting that the long time 
decay in the canonical dynamics was produced by acti- 
vation over low energy barriers. The scenario for this 
behavior is quite clear: after a long waiting time, the 
system diffuses further in a flat energy landscape sur- 
rounded by barriers of minimum height AE — 2. At 
time scales of the order of exp (2/T) it can further relax 
by thermal activation over these low barriers. The next 
level is at r ~ exp (4/T) « 2, 56 • 10 12 and it is clearly 
unreachable. The simplicity of the landscape seen by the 
system at this temperature explains the presence of the 
rapidly interrupted aging which is observed. 

In figure 4 we present the results of a simulation per- 
formed at a higher temperature (though still very low) 
T = 0.15T C . At this temperature in the time range of the 
simulation the system is able to probe activation over 
barriers of height AE = 2 and AE = 4. The waiting 
times are t w = 5000, 10000, 50000 and 200000. The three 
curves for the shortest waiting times can still be well col- 
lapsed for large values of t by the interrupted aging form 
(6), with an exponent /i « 0.3. However, the last curve 
(t w = 200000) clearly departs from this scaling, signaling 
a crossover to another regime. It is still possible to see 
an intermediate regime for the shortest t w 's. 



B. Full aging dynamics 

As the thermal energy is raised above the lowest lying 
levels the full ruggedness of the landscape should emerge 
and true aging dynamics should be restored. What hap- 
pens when the thermal energy is enough to turn discrete 
energy levels undetectable? In figure 5 we see the au- 
tocorrelations for a temperature T = 0.5T C for three 
waiting times t w = 10000, 50000, 100000. A fit to the 
quasi-equilibrium region gives x{T) = 0.1 and q = 0.62, 
indicating a very slow relaxation. This is due to the in- 
creasing complexity of the phase space visited at this 
temperature. The best scaling form obtained for the ag- 
ing regime is presented in figure 6 and corresponds to the 
following logarithmic form: 
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FIG. 4: Interrupted aging scaling for the asymptotic regime 
of the d=4 Edwards- Anderson model at T = 0.15T C with 
linear size L = 10 and four different waiting times t w = 
5000, 10000, 50000 and 200000. The fitting parameters of the 
stationary regime are A = 0.03 and x(T) = 0.54. In this 
regime the sub-aging exponent is /j, = 0.3. Note that the 
curve for the longest waiting time is beginning to depart from 
this regime. 



Note that there is still an intermediate regime which is 
not well described by this logarithmic scaling but the 
corresponding time window is too small to try a reason- 
able collapse there. On the other side, in the long time 
regime the logarithmic scaling (7) works remarkably well. 
It is worth noting that this form of the scaling function is 
compatible with dynamical ultrametricity (see 15 ). This 
form is slightly different from that expected by a droplet 
like scenario, which takes the form In (t + t w ) / In (t w ) and 
does not obey ultrametricity. This new evidence for a 
(weak) ultrametricity obtained directly from aging mea- 
surements is in agreement with recent results for the same 
model obtained from a quite different approach 16 . We 
have also done simulations at temperature T = 0.75T C 
(not shown) and the results are compatible with the log- 
arithmic scaling (7) although due to the strong noise at 
this rather high temperature the data is not so clean as 
at lower temperatures and the scaling of the curves is 
not so good. At high temperatures several time scales 
are mixed or superposed because of thermal noise and 
consequently it is difficult to obtain reliable scalings. 



C. d=6 

We have done a similar analysis of data obtained for 
the d — 6 Edwards- Anderson model. As this model is at 
the upper critical dimension we expected to see a behav- 
ior similar to the mean field or Sherrington-Kirkpatrick 
(SK) model 14 ' 17 ' 18 . The off equilibrium dynamics of the 
SK model is only poorly understood due to the complex- 
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FIG. 5: Two-time autocorrelation function for the ±J spin 
glass in d = 4 with linear size L = 10 after a quench to 
T = 0.5T C . The waiting times are t w = 10000, 50000 and 
100000. 



0.5625 



0.375 



0.25 




D = 4 T = 0.5T C 

^ = 10000 
L = 50000 
^ = 100000 



0.1 



log(t)/log(t w 



FIG. 6: Logarithmic aging scaling for the asymptotic time 
regime of the data shown in figure 5. Fit parameters are 
A = 0.3 and x(T) = 0.1. 



ity which emerges from the full replica symmerty break- 
ing which is the central characteristic of that model. As 
a consequence of the complex ultrametric organization of 
time scales no simple scaling form can be found for the 
aging dynamics which might be described by a superpo- 
sition of scaling regimes of the form 14 : 



and d = 4. The transition to an SK like aging scenario 
seems to be very slow as the connectivity of the system 
grows. 

All these simulations were carried out for systems with 
linear size L = 5. It is worth noting that the compu- 
tational time required for simulating these systems in- 
creases as L 6 , making it much more difficult to use larger 
values of L. The critical temperature of this model was 
estimated in 19 ' 20 to be T c ~ 3. Here again we start by 
considering the very low temperature case T = 0.07T C . 
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FIG. 7: Sub-aging scaling for two-time autocorrelation func- 
tions of the d=6 Edwards-Anderson ±J spin glass at T — 
0.07T C . The fit parameters are A = 0.027, x(T) = 0.55 and 
^ = 0.25. 
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FIG. 8: Sub-aging scaling for the two-time autocorrelation 
/o\ functions for the d=6 ±J spin glass at T = 0.15T C . Fit pa- 
V ' rameters are A = 0.04, x(T) = 0.29 and fi = 0.93. 



Nevertheless we have found a much simpler scenario 
which is in fact very similar to what is observed in d = 3 



The overall behavior of the two-time autocorrelation 
function is similar to that described for dimension d = 4. 
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It is possible to identify two different sub-aging regimes 
within the time window of our simulation, which covers 
up to t = 10 7 . Both in the intermediate and in the 
asymptotic regimes the functions C(t w + t,t w ) can be 
very well collapsed by an interrupted aging scenario of the 
form (6). In Figure 7 we show the best scaling obtained 
for the long time behavior of the autocorrelation after 
subtraction of the stationary part. The plots correspond 
to the three longest waiting times: t w — 2 15 , 2 17 and 
2 19 . The fitting in the quasi-equilibrium regime yielded 
x(T) = 0.55 and q ~ 0.96, similar to those found for d — 
4 at the same temperature. The two regimes described 
above can be clearly identified in this graph. For the 
last regime, the best data collapse yields an exponent 
fi w 0.25 and in the intermediate regime fi » 0.9 (this 
figure is not shown). 

In figure 8 we present the results of a similar simula- 
tion performed at T = 0.15T C . Here we note that the long 
times regime which had a smaller exponent (faster relax- 
ation) has now disappeared. Actually the three curves 
can be very well collapsed by a unique interrupted scal- 
ing form with fj, = 0.93. The main difference between this 
simulation and that performed for d = 4 (see figure 4) is 
the fact that now it is not seen the crossover described 
in d = 4 for the longest waiting time. 

Note also that as T increases the scaling form ap- 
proaches the simple aging scaling law /it — ► 1. 

At T = 0.5T C the scenario is not so clear as in d = 4. 
In figures 9 and 10 we present two different scalings. In 
9 we use the expression (7) and in 10 the same scaling 
form used in the interrupted scenario. Nevertheless, it 
is important to stress that in this last case the exponent 
that produces the best collapse is // w 1.2, indicating 
the possible existence of a super-aging regime. However, 
as can be concluded from these figures, it is difficult to 
decide which is the true scaling law. In summary, as 
mentioned at the beginning of this section the d = 6 
model presents an aging dynamics very similar to the 
d = 4 case at corresponding temperatures and it seems 
to be still far from the limit of the mean field system 
studied in 17 ' 18 . 



IV. DISCUSSION AND CONCLUSION 
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FIG. 9: Logarithmic scaling for the two-time autocorrelation 
functions of the d=6 Edwards- Anderson spin glass at T = 
0.5T C . Fit parameters are A = 0.34 and x(T) = 0.16 and 
q = 0.59. 
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FIG. 10: Super-aging scaling for the same data of figure (9). 
Fit parameters are A = 0.34 and x(T) = 0.16, q = 0.59 and 
/i=1.2. 



In this paper we have presented an extensive numerical 
study of the out of equilibrium dynamics of spin glasses 
with discrete couplings defined on hypercubic lattices in 4 
and 6 dimensions, for very low temperatures T — 0.07T C , 
T = 0.15T C and T = 0.5T C . 

The relaxation behavior of the model with ± J cou- 
plings is different from that with continuous gaussian 
couplings. While the observed aging dynamics of the 
gaussian model is described in a wide range of temper- 
atures and time scales by simple aging or weak inter- 
rupted aging scalings, the observed behavior of the ±J 
model seems to depart from that. At the lowest temper- 
atures studied , that is for T = 0.07T C and T = 0.15T C , 



the thermal energy is still not enough to permit activa- 
tion over several scales and only the lowest lying energy 
levels are probed. As a consequence the system diffuses 
on a flat energy landscape surrounded by low barriers of 
height 2 or at most 4. It relaxes in this simple landscape 
presenting a kind of interrupted aging with a rapid relax- 
ation to stationarity. In this sense these relaxations are 
practically independent of t w and true aging is restored 
only at higher temperatures. In order to see true aging at 
these very low temperatures we would have to wait for ex- 
tremely long times. At T = 0.5T C both in d = 4 and d = 6 
the full aging is restored, no more signs of interrupted 
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scalings are seen but a slower logarithmic scaling is ob- 
served. At this intermediate temperature range thermal 
energy is large compared with the low lying energy states 
which were important at the lower temperatures. Con- 
sequently the aging dynamics proceeds slowly but with 
the tendency of restoring ergodicity and equilibrium as t w 
grows. It is important to stress that for no one of the tem- 
peratures studied T = 0.07T C , T = 0.15T C , T = 0.5T C , 
and T = 0.75T C simple aging scenarios were found, in 
contrast to what seems to be the rule for the gaussian 
spin glass model. Insted the logarithmic aging observed 
at T = 0.5T C (and also at T — 0.75T C ) points to a simple 
hierarchy of time scales with dynamical ultrametricity as 
observed in 16 . On the other hand, even for the model in 
d = 6, which is at the upper critical dimension, the qual- 
itative form of the aging curves and the scaling functions 
are very different from those found in the SK model 17 
and in the Hopfield model 18 which are known to present 
full replica symmetry breaking and a consequently full 
hierarchical or ultrametric organization of time scales 14 . 
In this respect it is worth citing recent work by Yoshino, 
Hukushima and Takayama 12 where they presented an ex- 
tended version of the dynamical droplet theory. Two in- 
teresting new predicitions are the presence of a new dy- 
namical order parameter qo < Qea which should be ob- 
served in a particular lenght/time scaling regime and the 
distinction between the times at which time translation 



invariance and the Fluctuation-Dissipation Theorem are 
violated. The authors presented numerical results on the 
d = 4 ± J spin glass to support their findings but clearly 
more precise measurements are needed in order to test 
the theoretical predictions. In particular, scaling func- 
tions depend on the growth law of the coherence lenght 
L{t) which permits to see the crossover between critical 
to activated dynamics. Unfortunately this information is 
still not available for the temperatures and time scales 
reached in our simulations which we think are necessary 
in order to clearly see the separation of time regimes dur- 
ing aging dynamics. We expect new interesting results 
to come from these kind of analysis in the near future. 
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